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Abstract 

A search for direct pair production of the supersymmetric partner of the top quark, decaying 
via a scalar tau to a nearly massless gravitino, has been performed using 20 fb”^ of proton- 
proton collision data at = 8 TeV. The data were collected by the ATLAS experiment 
at the LHC in 2012. Top squark candidates are searched for in events with either two had- 
ronically decaying tau leptons, one hadronically decaying tau and one light lepton, or two 
light leptons. No significant excess over the Standard Model expectation is found. Exclusion 
limits at 95% confidence level are set as a function of the top squark and scalar tau masses. 
Depending on the scalar tau mass, ranging from the 87 GeV LEP limit to the top squark 
mass, lower limits between 490 GeV and 650 GeV are placed on the top squark mass within 
the model considered. 
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1 Introduction 

Additional partners of the top quark are ingredients in several models that address the hierarchy prob¬ 
lem [1^] of the Standard Model (SM). Supersymmetry (SUSY) [5-13] is one such model which natur¬ 
ally resolves the hierarchy problem with the introduction of supersymmetric partners of the known bosons 
and fermions. A supersymmetric partner of the top quark would stabilise the Higgs boson mass against 
quadratically divergent quantum corrections, provided that its mass is close to the electroweak symmetry 
breaking energy scale. This would make its discovery possible at the Large Hadron Collider (LHC) [14]. 
In a generic 7?-parity-conserving minimal supersymmetric extension of the SM (MSSM) [15-19], the 
scalar partners of right-handed and left-handed quarks, and q^, can mix, as can the scalar partners 
of charged leptons, and to form two squark or two slepton mass eigenstates, respectively. The 
lighter of the two top squark eigenstates is denoted t\ and is referred to as the scalar top in the following. 
Likewise, the lighter of the two scalar tau eigenstates is denoted f i and referred to herein as the scalar 
tau. 

In gauge-mediated supersymmetry breaking (GMSB) models [20-25], the spin-3/2 partner of the grav¬ 
iton, called the gravitino G, is assumed to be the lightest supersymmetric particle. Assuming that the 
mass scale of the messengers responsible for the supersymmetry breaking is of the order of 10 TeV, in 
order to minimise fine tuning [26], the scalar top should be lighter than about 400 GeV [27]. If the scalar 
tau is lighter than the scalar top, and the supersymmetric partners of the gauge and Higgs bosons (char- 
ginos and neutralinos) are heavier, the dominant decay mode of the t\ might be the three-body decay 
into bvrf\, where Vr is the tau neutrino, followed by the fi decay into a tau lepton and a gravitino. The 
other possible decay mode is the two-body decay into a top quark and a gravitino. The partial width of 
the two-body decay depends on the gravitino mass, while the partial width of the three-body decay via a 
virtual chargino depends on the chargino mass, as well as the chargino and scalar top mixing. For fixed 
scalar top and scalar tau masses either mode can dominate, and we focus in this paper on the signature 
resulting from the three-body decay. The two-body decay would give a signature very similar to that 
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of the decay into a top quark and a neutralino, which has been addressed in previous searches [28-34]. 
In the simplest gauge-mediated models, the predicted Higgs boson mass [35] is typically lower than the 
measured mass [36], especially if a light scalar top is also required. However, a variety of mechanisms 
exist [37-41] to raise the Higgs boson mass to make it compatible with the observed value. 

A lower limit of 87 GeV on the mass of the scalar tau has been set by the LEP experiments [42-46]. No 
limits have been published so far from hadron collider searches for the three-body decay of the scalar 
top into the scalar tau. Searches for scalar top pair production in proton-proton (pp) collisions, targeting 
the decay into charginos or neutralinos, have been performed by the ATLAS [28] and CMS [29-34] 
collaborations. Searches for scalar tops decaying into gravitinos, but not including the scalar tau in the 
decay chain, have been reported by the ATLAS [47] and CMS [48, 49] collaborations. 

This paper presents a dedicated search for pair production of scalar tops resulting in a final state with 
two tau leptons, two jets that contain a fj-hadron (fj-jets), and two very light gravitationally interacting 
particles. The decay topology of the signal process is shown in Lig. 1 ; the model considered is a simplified 
model in which all fhe supersymmefric particles ofher fhan fhe scalar lop and Ihe ones enfering ils decay 
chain are decoupled. In order fo maximise fhe sensilivily, fwo dislincl analyses have been performed 
based on fhe decay mode of fhe fau leplons in fhe final slafe: one analysis requires fwo hadronically 
decaying fau leplons (fhe hadron-hadron channel) and fhe ofher requires one hadronically decaying fau 
leplon and one fau decaying info an eleclron or muon, plus neulrinos (fhe leplon-hadron channel). In 
addilion, fhe resulls of fhe search reporled in Ref. [50], which is sensitive lo evenfs where bolh fau leplons 
decay leplonically (referred lo as fhe lepfon-lepfon channel), are reinlerpreled and limifs are sel on fhe 
scalar lop and scalar fau masses. 
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Figure 1: Diagram showing the decay topology of the signal process. 
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2 ATLAS detector 


ATLAS [51] is a multi-purpose particle physics experiment at the LHC. The ATLAS detector^ consists 
of an inner tracking detector surrounded by a superconducting solenoid, electromagnetic and hadronic 
calorimeters, and a muon spectrometer. The inner detector covers \r]\ < 2.5 and consists of a silicon 
pixel detector, a semiconductor microstrip detector, and a transition radiation tracker (TRT). The inner 
detector is surrounded by a thin superconducting solenoid providing a 2 T axial magnetic field, and allows 
for precision tracking of charged particles and vertex reconstruction. The calorimeter system covers the 
pseudorapidity range \t]\ < 4.9. In the region \r}\ < 3.2, high-granularity liquid-argon electromagnetic 
sampling calorimeters are used. A steel/scintillator-tile calorimeter provides energy measurements for 
hadrons within \r]\ < 1.7. The end-cap and forward regions, which cover the range 1.5 < I 77 I < 4.9, 
are instrumented with liquid-argon calorimeters for electromagnetic and hadronic particles. The muon 
spectrometer surrounds the calorimeters and consists of three large superconducting air-core toroid mag¬ 
nets, each with eight coils, a system of tracking chambers (covering I 77 I < 2.7) and fast trigger chambers 
(covering \j]\ < 2.4). 


3 Monte Carlo simulations and data samples 

A number of Monte Carlo (MC) simulated event samples are used to model the signal and describe the 
backgrounds. For the main background components, predictions are normalised to the data in control 
regions (CRs) and then extrapolated to the signal regions (SRs) using simulation. All MC samples util¬ 
ised in the analyses are processed using either the ATLAS detector simulation [52] based on GEANT4 [53] 
or a fast simulation based on a parameterisation of the performance of the ATLAS electromagnetic and 
hadronic calorimeters [54] and GEANT4 elsewhere. Additional pp interactions in the same (in-time) and 
nearby (out-of-time) bunch crossings, termed pile-up, are included in the simulation, and events are re¬ 
weighted so that the distribution of the number of pile-up collisions matches that in the data. 

The signal model considered is a supersymmetric model with the gravitino as the lightest supersymmetric 
particle. By construction, the scalar partner of the right-handed tau lepton and the lightest scalar top^ are 
the next-to-lightest and the next-to-next-to-lightest supersymmetric particles, respectively, and different 
signal models are simulated by varying their masses. Pair production of the scalar top is generated us¬ 
ing HERWIG-i-h 2.6.3 [55] with the parton distribution functions (PDF) set CTEQ6L1 [56]. The model 
requires that the scalar top decays to bvTf\ via a virtual chargino with 100 % branching ratio, while the 
f 1 decays, with a 100% branching ratio, into a tau lepton and a gravitino. Lifetimes are assumed to be 
small enough (below about 1 ps) that the detector response is unaffected by the decay distance of the 
supersymmetric particles from the primary vertex. 

Signal cross sections are calculated to next-to-leading order (NLO) in the strong coupling constant as, 
adding the resummation of soft gluon emission at next-to-leading-logarithmic accuracy (NLO-i-NLL) [57- 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis coinciding with the axis of the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the 
^-axis points upwards. Cylindrical coordinates (r, are used in the transverse plane, (j) being the azimuthal angle around the 
beam pipe. The pseudorapidity is defined in terms of the polar angle 6 as 77 = - In tan(6/2). 

^ The mixing matrix of the simulated samples is such that the lightest scalar top eigenstate is almost a pure partner of the 
right-handed top quark. 
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Table 1; Details about the MC generation of the background and signal samples. 


Process 

Generator 

Parton shower 

Cross-section 

normalisation 

PDF set 

Generator 

tune 

tt 

POWHEG-BOX 
r2129 [61, 62] 

PYTHIA 6.426 [63] 

NNLOh-NNLL 

[64-69] 

NLO CTIO 
[70] 

Perugia 
2011C [71] 

Single-top 

POWHEG-BOX 

PYTHIA 6.426 

NNLOh-NNLL 

CTEQ6L1 

Perugia 

iWt and s-channel) 

rl556 [61, 72, 73] 

[74] 

[56] 

2011C 

Single-top 

(f-channel) 

ACERMC 3.8 [75] 

PYTHIA 6.426 

NNLOh-NNLL 

[76] 

CTEQ6L1 

Perugia 

2011C 

tt + wiz 

MDGRAPH5 
1.3.28 [77] 

PYTHIA 6.426 

NLO [78] 

CTEQ6L1 

AUET2 

[79] 

IVIV, IVZ, zz 

SHERPA 1.4.1 [80] 

SHERPA 1.4.1 

NLO [81] 

NLO CTIO 

SHERPA 

default 

Zly*{-^ ee/pp)-Hjets 

ALPGEN 2.14 [82] 

HERWIG 6.520 [83] 

NNLO [84] 

CTEQ6L1 

AUET2 

Zly*{-^ TT)-Hjets 

SHERPA 1.4.1 

SHERPA 1.4.1 

NNLO [84] 

NLO CTIO 

SHERPA 

default 

W(-^ ^v)4-jets, 

£ — e,p,T 

SHERPA 1.4.1 

SHERPA 1.4.1 

NNLO [84] 

NLO CTIO 

SHERPA 

default 

ht\ 

HERWIG++ 2.6.3 
[55] 

HERWIG++ 2.6.3 

NLOh-NLL 

[57-59] 

CTEQ6L1 

UE-EE-3 

[85] 


59]. The nominal cross section and its uncertainty are taken from an envelope of cross-section predictions 
using different PDF sets and factorisation and renormalisation scales, as described in Ref. [60]. 

The programs used to generate signal and background events, as well as details of the cross-section 
calculation, PDF sets, and generator tunings, are reported in Table 1. 

The data sample used in this paper was recorded between March and December 2012, with the LHC 
operating at a centre-of-mass energy of = 8 TeV. The data are collected based on the decisions of a 
three-level trigger system [86]. Events are selected for the electron-hadron channel if they are accepted by 
a single-electron trigger, and for the muon-hadron channel if accepted by a single-muon trigger. For the 
hadron-hadron channel, a missing transverse momentum trigger is used. The trigger efficiency reaches 
its maximum value for leptons with a transverse momentum (pt) above 25 GeV in the lepton-hadron 
channels, and it exceeds 97% for a missing transverse momentum above 150 GeV in the hadron-hadron 
channel. After beam, detector and data-quality requirements, the integrated luminosity of the data samples 
is 20.3 fb“^ in the electron-hadron and muon-hadron channels, and 20.1 fb“^ [87] in the hadron-hadron 
channel. The difference in integrated luminosity is due to the additional data-quality requirements related 
to the trigger used in the hadron-hadron channel. 


4 Event reconstruction 

The reconstruction and selection of final-state objects used in this analysis are discussed below. 

Vertex candidates from pp interactions are reconstructed using tracks in the inner detector. To identify 
the hard-scattering vertex in the presence of pile-up, the vertex with the highest scalar sum of the squared 
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transverse momentum of the associated tracks, is defined as the primary vertex. The primary vertex 
is required to have at least five associated tracks with pT > 400 MeV. 

Jets are reconstructed from three-dimensional clusters of energy deposits in the calorimeters using the 
anti-kf jet clustering algorithm [ 88 ] using Fast Jet [89], with a radius parameter of /? = 0.4. The dif¬ 
ferences in the calorimeter response between electrons/photons and hadrons are taken into account by 
classifying each cluster as coming from a hadronic or an electromagnetic shower on the basis of its 
shape [90]. The energy of electromagnetic and hadronic clusters is then weighted with correction factors 
derived from MC simulations. The average expected contribution from pile-up, calculated as the product 
of the jet area and the median energy density of the event [91], is subtracted from the jet energy. A further 
energy and r] calibration based on MC simulations and data, relating the response of the calorimeter to 
the true simulated jet energy [92, 93], is then applied. The jets selected in the analysis are the jet candid¬ 
ates with pj > 20 GeV and \r]\ < 2.5. Events containing jets that are likely to have arisen from detector 
noise, beam background or cosmic rays, are removed using the procedures described in Ref. [92]. Events 
containing any jet failing to meet specific quality criteria described in Ref. [94] are also rejected. 

Among the jets satisfying the selection criteria above, b-]et candidates are identified by a neural-network- 
based algorithm, which utilises the impact parameters of tracks, secondary vertex reconstruction, and the 
topology of b- and c-hadron decays inside a jet [95, 96]. The efficiency for tagging f?-jets in a MC sample 
of tt events using this algorithm is 70% with rejection factors of 137 and 5 against light-quark or gluon 
jets, and c-quark jets, respectively. To compensate for differences between the f>-tagging efficiencies and 
mis-tag rates in data and MC simulation, correction factors derived using tt events are applied to jets in 
the simulation as described in Refs. [95, 96]. 

Electron candidates used to veto events with prompt leptons in the hadron-hadron channel search are 
required to have p^ > 10 GeV, I 77 I < 2.47 and to satisfy loose selection criteria on electromagnetic shower 
shape and track quality [97]. Their longitudinal and transverse impact parameters must be within 2 mm 
and 1 mm of the primary vertex, respectively. In the lepton-hadron channel, further selections are applied. 
Electrons are required to satisfy the tight quality criteria, to have pT > 25 GeV, and to be isolated within 
the tracking volume. The electron identification efficiencies are of about 95%, 91% and 80% for the 
loose, medium and tight working points respectively. The electron isolation requires that the scalar sum, 
EpT^ of the pt of inner detector tracks within a cone of size AR = '\/(A?7)2"T(A0)2 = 0.2 around the 
electron candidate, is less than 10% of the electron pj. The tracks included in the scalar sum must have 
Pt > 1 GeV, are matched to the primary vertex, and do not include the electron track. 

Muon candidates are reconstructed using inner detector tracks either combined with muon spectrometer 
tracks or matched to muon segments [98]. They are required to have px > 10 GeV and \r]\ < 2.4. Their 
longitudinal and transverse impact parameters must be within 1 mm and 0.2 mm of the primary vertex, 
respectively. These selections have an overall efficiency of about 99%. Muon candidates that pass these 
selections are referred to as loose muons and are used to veto events with prompt leptons in the hadron- 
hadron channel. The candidates with px > 25 GeV which fulfill the isolation requirement Epx <1.8 GeV, 
i.e. with at most one additional track with 1 < px < 1.8 GeV reconstructed within a cone of size AR = 0.2 
around the muon track, are referred to as tight muons. 

Event-level weights are applied to MC events to correct for differences between the lepton reconstruction 
and identification efficiencies measured in the simulation, and those measured in data. 

Hadronically decaying tau lepton (rhad) candidates are seeded by calorimeter jets with px > 10 GeV. An 
p- and px-dependent energy scale calibration is applied to correct for the detector response and subtract 
energy from pile-up interactions [99]. Tau lepton candidates are identified by using two boosted decision 
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Table 2: Sequence of the overlap removal algorithm. Here, { refers to electrons and muons. 


Condition 

Discarded objecf 

A/?(jel, elecfron) < 0.2 

jet 

ARfThad, i) < 0.2 

'^'had 

AR(jef,^) < 0.4 

€ 

AR(Thad,jet) < 0.2 

Jet 


tree (BDT) algorithms that separate them from jets and electrons [99]. Variables describing the shower 
shape in the calorimeters and information from the tracking system are used to separate the collimated Thad 
decay products from the generally broader jets resulting from quark and gluon hadronisation. Variables 
such as the number of tracks or the fraction of the total tau energy contained in a cone of size A/? = 0.1 
centred on the tau candidate provide strong discriminating power. To distinguish taus from electrons, 
the most discriminating characteristics are the transition radiation emitted by electrons in the TRT and 
the longer and wider shower generated by a hadronically decaying tau in the calorimeter compared with 
that produced by an electron. In addition to the two BDT selection criteria, a muon veto is applied. 
Hadronically decaying tau lepton candidates are required to have pj > 20 GeV, I 77 I < 2.47, and exactly 
one or three associated inner detector tracks (referred to as 1-prong and 3-prong candidates, respectively). 
The tau candidate is assigned an electric charge equal to the sum of the charges of the associated tracks, 
and this is required to be either -i-l or -1. Three working points {loose, medium, and tight) are used for 
each BDT. The hadron-hadron channel uses the tight identification working point for jet rejection and 
the medium identification working point for electron rejection, while the lepton-hadron channel uses the 
medium working point for both. The loose working point has been used to cross-check the background 
modelling. For the jet-veto BDT, the working points correspond to a signal efficiency of 70%, 60% and 
40% for 1-prong Thad, and 65%, 55% and 35% for 3-prong Thad, respectively. The electron-veto BDT 
working points have a signal efficiency of 95%, 85% and 75%, respectively. Efficiency scale factors 
are used to account for the mis-modelling of BDT input variables in the simulation. They are extracted 
by comparing efficiencies in data and simulation in a Z ^ tt selection, using a tag-and-probe method 
described in Ref. [99]. 

As a given final-state particle can be simultaneously reconstructed as (for example) an electron, a jet and 
a hadronically decaying tau lepton, an algorithm is used to resolve such ambiguities. Electrons satisfying 
the medium quality criteria, muons satisfying the criteria described above except that on isolation, jets 
and hadronically decaying tau candidates satisfying the selection criteria given above are considered by 
the algorithm. If two objects are close together in AR, one of them is discarded according to the sequence 
specified in Table 2. Elecfrons and muons close fo jefs, which are likely fo originafe from fhe decay of 
heavy-flavour hadrons, are removed from fhe lisf of lepfons used in fhe analysis. 

The missing fransverse momenfum vecfor p™*®, whose magnifude is referred fo as is calculafed 

as fhe negafive vecfor sum of fhe fransverse momenfa of all reconsfrucfed elecfrons, jefs and muons, and 
calorimefer energy clusfers nof associafed wifh any objecfs. Eor fhe compufafion, hadronically de¬ 
caying fans are freafed as jefs. Clusfers associafed wifh elecfrons wifh pT > 10 GeV, and fhose associated 
wifh jefs wifh pT > 20 GeV are calibrated wifh fhe elecfron and jef clusfer calibrations, respectively. Eor 
jefs, fhe calibration includes fhe pile-up correction described earlier while the. jet vertex fraction (JVE) 
requiremenf is nof imposed. The JVE variable is fhe rafio of fhe sum of fhe fransverse momenfum of fhe 
fracks associafed wifh fhe jef and originating from fhe selecfed primary vertex fo fhe fofal pT sum of all 
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tracks matched with the jet. This requirement rejects jets originating from pile-up. Clusters of energy 
deposits in calorimeter cells with \r}\ < 2.5 not associated with these objects are calibrated using both 
calorimeter and tracker information [100]. 


5 Event selection and background estimate 

5.1 Hadron-hadron channel 

For the hadron-hadron channel search, events in the signal region are required to have exactly two oppos¬ 
itely charged hadronically decaying taus satisfying the tight identification criteria, no electrons or muons, 
and at least two jets with a JVF larger than 0.5 or pj > 50 GeV. One of the jets must be fj-tagged. The 
leading jet must also satisfy pj > 40 GeV. 

The missing transverse momentum must be larger than 150 GeV. The A0 separation between each of the 
two leading jets and the direction of the missing transverse momentum must be greater than 0.5 radian, to 
suppress events where large F"™®® arises from mis-measurement of jet energies. Beyond these preselection 
requirements, additional selections are made using transverse masses and derived variables, as explained 
below. These selections have been determined using MC signal and background samples to maximise the 
expected significance of fhe signal. 

The fransverse mass associafed wifh fwo final-slafe objecfs a and b is defined as 

mT(a,b)= yjml + m^ + 2 , (1) 

where m, E^ and pj are fhe objecf mass, fransverse energy and fransverse momenfum vecfor, respecfively. 
Objecfs enfering fhe mj calculafion are always assumed fo be massless, unless fhe fransverse mass is used 
as part of a derived variable in fhe leplon-hadron channel (see Seel. 5.2). 

The stransverse mass {mj 2 ) [101, 102] is compufed as 


mji (a, b) 



( 2 ) 


where and q!j, are veefors salisfying q^ + q!j. - p™®®, and fhe minimum is faken over all fhe possible 
choices of q“ and q!j,. 

The selecfion criferia fhaf define fhe signal region for fhe hadron-hadron channel (SRHH) rely on fhe 
following variables: 


• mT 2 (Thad> Thad) is defined using fhe momenla of fhe hadronically decaying fans and fhe missing 
fransverse momenfum, which is assumed lo resull from fwo invisible massless particles. The 
mT 7 (T hari , Thad) variable is bounded from above by fhe W boson mass for evenfs where fhe fwo had¬ 
ronically decaying fans originate from fhe decay of fwo W bosons and all fhe missing fransverse 
momenfum is carried by fhe neulrinos from fhe W bosons decay, as is fhe case for fhe dominanl 
background (tt). 




• wJ™™(Thad, "rhad) is defined as the sum of the transverse mass of each Thad candidate and the missing 
transverse momentum 


(Thad, Thad) ^ (rhadl, P™") + (T'had2, P™'') (3) 

The m™™(Thad, Thad) distribution is expected to reach higher values for the signal due to a larger 
number of invisible final-state particles than for the SM background processes. 

For the SRHH signal region, the stransverse mass /nxiCThadi,Thad 2 ) is required to be larger than 50 GeV 
while the ni™™(Thadi, Thad 2 ) variable is required to be larger than 160 GeV. The signal selection efficiency, 
defined as the number of signal events that pass the full selection over the total number of generated 
events, is only weakly dependent on the scalar tau mass, while it increases from 0.02% to 0.7% as the 
scalar top mass increases from 150 GeV to 700 GeV, for a scalar tau mass of 87 GeV. The distributions 
of /nT 2 (Thadi, Thad 2 ) and m™’"(Thadi, Thad 2 ) are illustrated in Fig. 2 after the preselection. 

The background processes populating the SRHH selection are grouped into three categories. The first 
contains events with two real, hadronically decaying taus {true taus). It consists mainly of tt events, 
with smaller contributions from single-top-quark, Z-rjets, diboson {WW, fkZ, ZZ) and tt + V production, 
where V = fT, Z. This set of backgrounds is estimated from simulation. The remaining backgrounds 
contain events where at least one tau candidate is an electron or a jet that passes the tau identification 
criteria (fake taus). The second category, which contains events with only one fake Thad, is composed of 
tt, single-top-quark and VT-i-jets events. The third and smaller category corresponds to processes with two 
fake taus. It is mostly composed of tt, Z(—> yy)-i-jets, and single-top-quark events, which are all estimated 
from simulation. It has been verified that these backgrounds are well modelled: in kinematic selections 
where ti with true taus is expected to be the dominant process, the ratio of data over the MC prediction 
is compatible with one within systematics uncertainties. The contribution from multi-jet events, where 
both tau candidates are fakes, is estimated from data using the jet smearing method described later in this 
section. 

The single-fake Thad backgrounds from top quark (tt and single-top) and VT-i-jets events are estimated 
using MC simulations scaled to the observed number of data events in two dedicated control regions 
(CRHHTop and CRHHWjets). These control regions require a single-muon trigger, one Thad satisfying the 
tight quality criteria, and one muon with px > 25 GeV that satisfies the tight quality criteria. The niT 2 and 
variables are then calculated using the tau and muon momenta, considering the invisible particles 
as massless. One muon and one Thad are required in the control regions rather than two hadronically 
decaying taus in order to minimise signal contamination. Upper bounds are set on the ^72 and 
variables, which make the contamination from the lepton-hadron signal negligible. Table 3 details the 
selections defining the two control regions and the signal region. The contributions to the background 
from the double-fake Thad sources are smaller than 4.5% and therefore they are estimated using simulation 
without normalising to data in a control region. 

A simultaneous likelihood fit is performed to determine the normalisation factors of the single-fake Thad 
backgrounds, with the number of data events in each CR as constraint, and the systematic uncertainties 
described in Sect. 6 included as nuisance parameters. The fit is used to predict the number of background 
events in the CRs and the SR. The background modelling is verified using two validation regions (VRs) 
by comparing the observed number of events in each VR with the number derived from the fit. The 
single-fake Thad backgrounds from top quark and IT-i-jets events each have a validation region, labelled 
VRHHTop and VRHHWjets. Like the control regions, they are defined using a muon and tau to avoid 
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Figure 2: Top; distribution of the stransverse mass constructed from the two Thad, 'WT 2 ('!‘had, Thad), for events passing 
the hadron-hadron preselection requirements. Bottom: distribution of the sum of the transverse mass of each Thad 
candidate and the missing transverse momentum, /«j“™(Thad,Thad), for events passing the hadron-hadron preselec¬ 
tion requirements. The contributions from all SM backgrounds are shown as a histogram stack; the bands represent 
the total uncertainty. The distributions expected for two signal models are also shown. 
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signal contamination, and the selections are summarised in Table 3. The validation regions are designed 
to be kinematically close to the signal region without overlapping with the control or signal regions. 
The composition of the control and validation regions after the fit is shown in Fig. 3. The observed and 
expected background yields in the VRs are in good agreement, with 50 observed events in VRHHWjets 
(48.5 + 6.9 expected) and 31 observed events in VRHHTop (29.0 ±4.1 expected). It has also been verified 
that a normalisation factor for the top quark background with two real Thad would be compatible with one 
within uncertainties. 

The multi-jet background is estimated from data using the jet smearing method described in Ref. [103]. 
A set of single-jet triggers is used to select a sample of events with at least two jets (of which at least 
one is required to be a b-]e.t), and two Thad candidates. These events are required to have a low 
significance^, fo refain topologies where jets and tau candidates are well-balanced in the transverse plane 
and suppress processes with genuine F™®®. The energy of jets and tau candidates is then smeared within 
the resolution of the calorimeter, in order to simulate F™®® arising from mis-measurements. To minimise 
the statistical uncertainty, no identification criteria are applied to Thad candidates beyond the 1,3-track 
requirement, and a fake rate is used at a later stage to account for the tau identification efficiency. The 
pseudo-datasef obtained after smearing serves as a template for the multi-jet background. Its normalisa¬ 
tion is derived in a multi-jet-enriched CR, labelled CRHHQCD in Table 3. To estimate the background 
yield in the signal region, all SRHH requirements except the tau identification are applied to the norm¬ 
alised background template. A weight is then applied to each event according to the probability for a jet 
reconstructed as a tau to satisfy the tight tau identification criteria. This fake rate is measured in data using 
events which fire a single-jef frigger, with at least two jets and a hadronically decaying tau candidate. It 
is found to be of the order of 1% for 1-prong tau candidates and between 0.02% and 0.4% (with a strong 
Pt dependence) for 3-prong tau candidates. The number of multi-jet events in the SR is estimated to be 
0.0043 ± 0.0007 (stat) Ig oqqs therefore neglected. 


Table 3: Definition of the signal region (SRHH) for the hadron-hadron analysis. The selections of the associated 
control regions for tt and single-top-quark (CRHHTop) and lT±jets (CRHHWjets) events with one fake hadronically 
decaying tau, as well as the validation regions (VRHHTop and VRHHWjets), are also shown. The ( entering the 
mj 2 and m™™ variables is either a Thad (SR) or a muon (CRs and VRs). 


Region 


N, 

fVjet 

^ fc’-jet 

^miss 

AcPUi.2,pT^) 

tMT2(Thad, i) 

mr(^had,Q 

SRHH 

2 

0 

>2 

> 1 

> 150 GeV 

>0.5 

> 50 GeV 

> 160 GeV 

CRHHTop 

1 

1 

>2 

> 1 

> 100 GeV 

>0.5 

- 

[70,120] GeV 

CRHHWjets 

1 

1 

>2 

0 

> 100 GeV 

>0.5 

< 40 GeV 

[80,120] GeV 

VRHHTop 

1 

1 

>2 

> 1 

> 120 GeV 

>0.5 

< 40 GeV 

[120,140] GeV 

VRHHWjets 

1 

1 

>2 

0 

> 120 GeV 

>0.5 

< 40 GeV 

[120,150] GeV 

CRHHQCD 

> 2“ 

0 

> 2 

> 1 

> 150 GeV 

< 0.5* 

- 

- 


“ For the multi-jet control region (CRHHQCD), no identification criteria are applied to tau leptons. 
* The A0 requirement only applies to the sub-leading jet p- 


^ The significance is defined as 


2 £t + 

jets 


Yj Ej where soft terms correspond to clusters of energy deposits in 

soft terms 


the calorimeter which are not associated with any reconstructed object. 
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Figure 3: Background yields and composition after the fit in the two CRs and the two VRs of the hadron-hadron 
channel analysis. Combined statistical and systematic uncertainties are shown as shaded bands. The observed 
number of events and the total (constrained) background are the same by construction in the CRs. 


5.2 Lepton-hadron channel 

The search in the lepton-hadron channel requires exactly one hadronically decaying tau, exactly one 
isolated electron or muon with pT > 25 GeV, and no further isolated electrons or muons with pj > 
10 GeV. The hadronically decaying tau and the lepton are required to have opposite electric charge. Each 
event must also contain at least two jets, where at least one of the two jets must have pj > 50 GeV, and 
at least one of the two must be b-tagged. 

After this common preselection, two different signal regions are defined fo target signal models with a 
scalar top mass large or small in comparison to the top-quark mass. These are referred to as the low-mass 
(SRLM) and high-mass (SRHM) selections in the following, and they have been optimised with respect 
to the expected significance of the signal. The selections for the two signal regions are summarised in 
Tables 4 and 5. The low-mass selection requires a second ^-jet. Three niT 2 variables are employed in the 
selections, with different choices of the two visible four-momenta used in the calculation from Eq. (2): 

• mj 2 i^, Thad) uses the momenta of the light lepton and the hadronically decaying tau. The missing 
transverse momentum is assumed to result from two invisible massless particles. The nrT 2 (f’, Thad) 
variable is bounded from above by the W boson mass for events where the light lepton, the hadron¬ 
ically decaying tau and the missing transverse momentum originate from the decay of a pair of W 
bosons, which is the case for most of the background (tt and Wt). The high-mass selection requires 
this variable to be large, because its distribution for signal models with heavy scalar taus and scalar 
tops peaks at higher values than for the top-quark-dominated SM background. 

• mi 2 {b{, ^Thad) is calculated using the two jets with the highest fj-tagging weight. One of them is 
paired with the light lepton and the other with the Thad- The four-momentum vectors of the two 
resulting particle pairs are then used in the mT 2 algorithm. The missing transverse momentum is 
assumed to be carried by two invisible massless particles. Eor tt events where the jet and the lepton 
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belong to the decay of the same top quark, this variable is bounded from above by the top-quark 
mass. Similarly, for signal events, the upper bound on this variable is the scalar top mass. A 
maximum-value cut is therefore used in the low-mass selection and a minimum-value cut in the 
high-mass selection. The calculation of the variable requires the resolution of a two-fold ambiguity 
in the pairing of the jets and the leptons. Only the pairings for which m{M) and m(Z7Thad) are both 
smaller than nit are considered.^ If exactly one pairing satisfies the condition, that pairing is used 
in the mT 2 calculation. If both pairings satisfy the condition, mT 2 is calculated for both pairings 
and the smaller value is taken. If no pairing satisfies fhe condifion, fhe evenf is considered fo have 
passed fhe mT 2 ib£, fjThad) selecfion for fhe high-mass signal region and fo have failed if for fhe 
low-mass signal region. 

• mj 2 {b£, b) is only used for fhe low-mass selection. The sysfem of one of fhe fj-jels and fhe lighf 
lepfon is considered as fhe firsf visible four-momenlum. Only pairings for which m{b€) < nit are 
considered. If neifher pairing salisfies fhe condition, fhe even! is discarded, while if bofh pairings 
do, fhe pairing which yields fhe smaller value of mj 2 {b£, b) is used. The invisible parficle associafed 
wifh fhis sysfem is assumed fo be massless. The ofher b-]e.t is fhe second visible sysfem used in 
fhe /7 It 2 calculation, and fhe mass of fhe associafed invisible parficle is sef fo fhe W boson mass, as 
fhe algorifhm fargefs tt evenfs where one lepfon from a W boson decay is nof defecfed or identified. 
For fhe dominanf fop-quark background, fhe mT' 2 {b£, b) variable is bounded from above by fhe fop- 
quark mass. This variable has a sofler disfribution for low-mass signal evenfs than the background, 
and a maximum-value cut of 100 GeV is applied. 

The distributions for mT' 2 {b£, ^Thad) and mj 2 i£, rhad) are illustrated in Fig. 4 after the preselection, showing 
the separation between two signal models and the SM background. The mi 2 {b£, ^Thad) variable is used 
to distinguish the scalar top signal from the dominant top-quark backgrounds for both the low-mass and 
high-mass selections. 

Another variable used in the selections is the ratio of the scalar sum of the transverse momenta of the 
two leading jets (Ht) to the effective mass, /Weff = -i- Hj + where and are the 

transverse momenta of the lepton and the hadronically decaying tau, respectively. This ratio, Hjlmeg, 
tends to be smaller for signal events because of the high number of invisible particles in the final sfafe, and 
if is required fo be less fhan 0.5. The high-mass selection also requires fhe missing transverse momentum 
to be larger than 150 GeV and meff to be larger than 400 GeV because the decay products of a high-mass 
scalar top would have large momenta. The low-mass selection requires {pj + Ini,t« > 0.2 because 

the difference befween fhe masses of fhe scalar lop and scalar fau is relafively small in comparison fo fhe 
difference befween fhe masses of fhe lop quark and fhe W boson. Finally, fhe ni'i{£, p™^^) variable is used 
fo distinguish evenfs wifh real fau leptons from evenfs wifh fake fau leptons in the dominant top-quark 
background, and to distinguish multi-jet events from IT-i-jets events. The definitions of the low-mass and 
high-mass SRs are summarised in Tables 4 and 5, respectively. 

The signal selection efficiency of the low-mass selection is between 0.008% and 0.01% for the models 
with a scalar top mass between 150 GeV and 200 GeV, which is the target of this selection. The signal 
efficiency of the high-mass selection increases with the scalar top mass. For a fixed scalar top mass, if 
increases wifh fhe scalar fau mass as fhe mj 2 i£, Thad) selecfion becomes more efficienl, up fo fhe region 
wifh m(ti) - m{f) < 50 GeV where fhe f^-jets become too soft to be efficiently detected. Outside this 

For top-quark pair production events where the lepton and the jet belong to the decay of the same top quark, the invariant 
mass has an upper bound at — m^, approximately 152 GeV. The algorithm tries to select pairs that satisfy this condition, 
loosened to account for the detector resolution. 


13 






Figure 4; Top; distribution of the stransverse mass constructed from the h-jet plus lepton and h-jet plus Thad, 
mj 2 {b{, hrhad), for events passing the lepton-hadron preselection requirements with the additional requirement of a 
second h-tagged jet. Bottom; distribution of the stransverse mass constructed from the momenta of the light lepton 
and the hadronically decaying tau, mjiit, Thad), for events passing the lepton-hadron preselection requirements. 
The contributions from all SM backgrounds are shown as a histogram stack; the bands represent the total uncer¬ 
tainty. The overflow bin in the mx 2 (hf’, hzhad) plot is hlled with the events that have no (M, hrhad) pairing satisfying 
m{b() < m, and m(hrhad) < ni,. The distributions expected for two signal models are also shown. 
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region, which is better targeted by the lepton-lepton channel, the efficiency of the high-mass selection 
varies between 0.0007% and 1% for a scalar top mass between 200 GeV and 700 GeV. 


Table 4; Definition of the signal region SRLM used in the low-mass lepton-hadron analysis. The selections 
of the associated control regions for top-quark events with true taus (CRTtLM), top-quark events with fake taus 
(CRTfLM), and VT-i-jets (CRWLM), and of the validation region (VRTLM) are also given. 


Region 

^ b-jet 



mT2ib£, b) 

mi2{bi, (>Thad) 


W^eff 

SRLM 

> 2 

<0.5 

> 0.2 

< 100 GeV 

< 60 GeV 

- 

- 

CRTtLM 

> 2 

- 

> 0.2 

< 100 GeV 

110- 160 GeV 

> 100 GeV 

- 

CRTfLM 

> 2 

- 

> 0.2 

< 100 GeV 

110- 160 GeV 

< 100 GeV 

- 

CRWLM 

0 

<0.5 

> 0.2 

- 

- 

> 40 GeV 

< 400 GeV 

VRTLM 

> 2 

>0.5 

> 0.2 

< 100 GeV 

60- 110 GeV 

- 

- 


Table 5: Definition of the signal region SRHM used in the high-mass lepton-hadron analysis. The selections 
of the associated control regions for top-quark events with ttue taus (CRTtHM), top-quark events with fake taus 
(CRTfHM), and VT-i-jets (CRWHM), and of the validation region (VRTHM) are also given. 


Region 

N b-jet 

^miss 


^x/wJeff 

mj2{bi, bThad) 

mT2i{, Thad) 


SRHM 

> 1 

> 150GeV 

> 400 GeV 

<0.5 

> 180GeV 

> 120 GeV 

- 

CRTtHM 

> 1 

> 150GeV 

> 400 GeV 

<0.5 

> 180GeV 

20-80 GeV 

> 120 GeV 

CRTfHM 

> 1 

> 150GeV 

> 400 GeV 

<0.5 

> 180GeV 

20-80 GeV 

< 120 GeV 

CRWHM 

0 

> 150GeV 

> 400 GeV 

<0.5 

- 

20-80 GeV 

40-100 GeV 

VRHM 

> 1 

< 150GeV 

> 400 GeV 

<0.5 

> 180GeV 

> 80 GeV 

- 


In the lepton-hadron channel, the ratio of real to fake hadronically decaying tan events depends on the 
background process. In IT-rjets events, the light lepton is always a real lepton from the W decay, due to 
the high reconstruction efficiency and purity of final-state electrons and muons, while the Thad is faked by 
a recoiling hadronic object. In tt and Wt events, the light lepton originates from the decay of one of the W 
bosons while the hadronically decaying tau candidate can be either a real or a fake tau. These processes 
(VT-tjets, tt, and Wt) are the main background sources and are estimated by MC simulation scaled to the 
observed data in three CRs for each SR. The CRs are enriched in either IT-i-jets, top-quark events with 
true hadronically decaying taus, or top-quark events with fake hadronically decaying taus (where the top- 
quark events include both single and pair production), and are used to derive normalisation factors for 
these three categories of background. For the low-mass selection SRLM, the true- and fake-tau top-quark 
backgrounds are controlled by CRTtLM and CRTfLM, while CRWLM controls the W-i-jets background. 
For the high-mass selection SRHM, the three control regions CRTtHM, CRTfHM and CRWHM are used 
to normalise the true- and fake-tau top-quark backgrounds and the VF-i-jets background. The CRs are 
defined in Table 4 for the low-mass selection and in Table 5 for the high-mass selection. The minor 
contribution from other background processes is estimated from simulation. 

A simultaneous likelihood fit is performed to obtain the three normalisation factors for each SR, using the 
observed number of data events in each CR as constraints, and with the systematic uncertainty sources 
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(described in Sect. 6) treated as nuisance parameters. The fit is used to predict the number of background 
events in the CRs and the SR. The validity of the background modelling is verified by using a validation 
region for each SR and comparing the observed number of events with the prediction from the fit. For the 
low-mass selection, the validation region VRLM is defined in Table 4, while the validation region VRHM 
is defined in Table 5 for the high-mass selection. 

The background composition and the observed number of events in each CR as well as in the VR and SR 
are shown in Fig. 5 for the low-mass selection and in Fig. 6 for the high-mass selection. The observed and 
expected background yields in the VRs are in good agreement, with 386 observed events for the low-mass 
selection (351 ± 84 expected) and 17 observed events in the high-mass selection (22 ± 5 expected). The 
expected background yields and observed number of events in the SRs are reported in Sect. 7. 
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Figure 5: Background yields and composition after the fit for the three CRs and the VR in the lepton-hadron channel 
low-mass selection. Combined statistical and systematic uncertainties are shown as shaded bands. The observed 
number of events and the total (constrained) background are the same by constmction in the CRs. 

The background estimate with fake hadronically decaying taus (either from top-quark or VT-i-jets events) 
is validated using an alternative method. The observed rate of events with a light lepton and a Thad with 
the same electric charge is scaled by the expected ratio of opposite-sign (OS) to same-sign (SS) events for 
the fake Thad backgrounds, which is estimated from MC simulation. Too few SS events are observed for 
the SRHM selection to make a meaningful prediction, so the method is only viable for the looser SRLM 
selection, for which it predicts 12 ± 6 events with fake hadronically decaying taus, in agreement within 
uncertainties with the sum of IT-i-jets and top-quark events with fake hadronically decaying taus obtained 
from the fit, which is 12 ± 5 events. 


6 Systematic uncertainties 


Various sources of systematic uncertainty affecting the predicted background yields in the signal regions 
are considered. The uncertainties are either computed directly in the SR when backgrounds are estimated 
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Figure 6; Background yields and composition after the fit for the three CRs and the VR in the lepton-hadron channel 
high-mass selection. Combined statistical and systematic uncertainties are shown as shaded bands. The observed 
number of events and the total (constrained) background are the same by construction in the CRs. 


from simulation, or propagated through the fit for backgrounds that are normalised in CRs. 

The dominant detector-related systematic uncertainties considered in these analyses are the jet energy 
scale and resolution [92], the Thad energy scale and BDT identification efficiency [99], and fhe ft-tagging 
efficiency [95, 96]. The energy scale and resolution of clusters in the calorimeter not associated with 
electrons, muons or jets, which affect the missing transverse momentum calculation, are also a source 
of systematic uncertainty. In all cases, the difference in the predicted background or signal between the 
nominal MC simulation and that obtained after applying each systematic variation is used to determine the 
systematic uncertainty on the background or signal estimate. Parts of the systematic uncertainties cancel 
when a background is estimated from a control region, but they do not cancel for processes normalised 
to their theoretical cross section. The remaining detector-related systematic uncertainties, such as those 
on lepton reconstruction efficiency and on the modelling of the trigger, are of the order of a few percent. 
A 2.8% uncertainty on the luminosity determination was measured using techniques similar to that of 
Ref. [87], and it is included for the normalisation of all signal and background MC samples. The signal 
uncertainties are between 10% and 15% for models close to the observed exclusion contour. 

Various theoretical uncertainties are considered for the modelling of the major SM backgrounds. In the 
case of top-quark contributions, the predictions of POWHEG-BOX are compared with those of MC@NL0-4. ©6 
to estimate the uncertainty due to the choice of generator. The difference in the yields obtained from 
POWHEG-BOX interfaced to PYTHIA and POWHEG-BOX interfaced to HERWIG is taken as the systematic 
uncertainty due to parton shower modelling, and the predictions of dedicated ACERMC-3.8 samples gen¬ 
erated with different tuning parameters are compared to give the uncertainty related to the modelling of 
initial- and final-state radiation (ISR/FSR). At NLO, contributions with an additional bottom quark in 
the final state lead to ambiguities in the distinction between the Wt process {gb Wtb) and top-quark 
pair production. All the Wt samples, generated using MC@NL0-4.06 and POWHEG-BOX, use the diagram 
removal scheme [104] to model this interference. The ACERMC-3.8 event generator is used to simulate 
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the WWb and WWbb final states at leading order (which include both the tt and Wt single-top-quark 
processes); the predictions of these ACERMC-3.8 samples are then compared to those of the nominal MC 
samples in order to assess the uncertainty on the background estimate from this interference. The uncer¬ 
tainties on VT-i-jets and Z-i-jets production are evaluated by studying the predictions of ALPGEN-2.14 with 
various choices of the renormalisation and factorisation scales. 

The impact of systematic uncertainties on the total background estimate in the different SRs is shown in 
Table 6. The table quotes, for each SR, the relative background uncertainty attributed to each source. 


Table 6: Summary of background estimates and the associated total uncertainties.The size of each systematic 
uncertainty is quoted as a relative uncertainty on the total background. A dash indicates a negligible contribution 
to the uncertainty. The individual uncertainties can be correlated, and thus do not necessarily sum in quadrature to 
the total relative uncertainty. 



SRHH 

SRLM 

SRHM 

Background events 

3.1 ± 1.2 

22.1 ±4.7 

2.1 ± 1.5 

Uncertainty Breakdown [%]: 

Jet energy scale and resolution 

17 

13 

2 

Tau energy scale 

9 

4 

3 

Cluster energy scale and resolution 

1 

2 

4 

f7-tagging 

2 

4 

2 

Top-quark theory uncertainty 

37 

11 

64 

VF-i-jets theory and normalisation 

- 

1 

19 

Simulation statistics 

20 

6 

21 

Top normalisation 

18 

6 

20 


Signal cross sections are calculated at NLO-tNLL with a total associated uncertainty between 14% and 
16% for scalar top masses between 150 GeV and 560 GeV. 


7 Results and interpretation 

The numbers of events observed in the hadron-hadron SR and in the two lepton-hadron SRs are reported 
in Table 7, along with the background yields before and after the background-only likelihood fit. In 
both the results and interpretation tables (Tables 7 and 8) the quoted uncertainties include all the sources 
of statistical and systematic uncertainty. Good agreement is seen between the observed yields and the 
background estimates. 

Figure 7 shows the distributions of m™™(Thad, "Thad) and niT 2 (Thad 5 ^had) for the hadron-hadron channel, 
for events satisfying all the SR criteria except that on the variable being reported in the figure. Figure 8 
shows mj 2 {b£, f^Thad) for the lepton-hadron low-mass selection and Thad) for the lepton-hadron 

high-mass selection for events satisfying all the corresponding SR criteria except those on the variable 
displayed in the figure. 

Upper limits at 95% confidence level (CL) on the number of beyond-the-SM (BSM) events for each SR 
are derived using the HistFitter program [105], with the CL^ likelihood ratio prescription as described 
in Ref. [106]. The limits are calculated for each SR separately, with the observed number of events. 
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Figure 7: Top: distribution of ?M^™(Thad, fhad) for the events passing all the hadron-hadron signal region require¬ 
ments, except that on m^™(Thad, Thad)- Bottom: distribution of mxaCrhad, ^had) for the events passing all the hadron- 
hadron signal region requirements, except that on mT 2 (rhad, Thad)- The contributions from all SM backgrounds are 
shown as a histogram stack; the bands represent the total uncertainty. The background yields have been rescaled by 
the post-fit normalisation factors. The arrows mark the cut values used to define the SRs. The distributions expected 
for two signal models are also shown. 
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Figure 8: Top; distribution of mj 2 (bi, ^Thad) for events passing all the lepton-hadron LM signal region requirements, 
except that on inj 2 ib£, bTh-^d). Bottom: distribution of Thad) for events passing all the lepton-hadron HM 

signal region requirements, except that on mj 2 {£, Thad)- The contributions from all SM backgrounds are shown as a 
histogram stack; the bands represent the total uncertainty. The background yields have been rescaled by the post-fit 
normalisation factors. The arrows mark the cut values used to dehne the SRs. The overflow bin in the mj 2 {b{, hrhad) 
plot is filled with the events that have for both pairings of m(bi) and m{bT\^^) at least one invariant mass larger than 
m,. The distributions expected for two signal models are also shown. 
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Table 7: Observed number of events and background fit results for the hadron-hadron SR and the two lepton- 
hadron SRs. Combined statistical and systematic uncertainties are given. The uncertainties between the different 
background components can be correlated, so they do not necessarily sum to the total background uncertainty. A 
dash indicates a negligible background contribution. The nominal expectations from MC simulation are given for 
comparison in the lower part of the table. 


Channel 

SRHH 

SRLM 

SRHM 

Observed events 

3 

20 

3 

Total (constrained) background events 

3.1 + 1.2 

22.1+4.7 

2.1 ± 1.5 

Top with only true tau(s) 

2.0+ 1.1 

8.2 ±3.9 


Top with at least one fake tau 

0.9 ±0.5 

9.8 ±4.5 


IT-njets 

0 

^•^^- 0.01 

2.2 ± 0.6 

0.4 ± 0.4 

Z/r*+jets 

0 04 

- 0.04 

1.9+ 1.1 

- 

tt+V 

0.04 +0.02 

- 

0.3 ±0.1 

Diboson 

0.14 + 0.02 

- 

- 

Expected background events before the fit 

3.7 

25.8 

2.2 

Top with only true tau(s) 

2.0 

11.5 

0.18 

Top with at least one fake tau 

1.4 

10.1 

1.1 

VT-njets 

0.01 

2.4 

0.65 

Z/r*+jets 

0.04 

1.9 

- 

tt+V 

0.04 

- 

0.27 

Diboson 

0.14 

- 

- 


Table 8; Left to right; Total constrained background yields, number of observed events, 95% CL observed (expec¬ 
ted) upper limits on the number of BSM events, , ,, and the visible cross section, (^e(T)^l , ,. 

' obs.(exp.)’ ’ ' 'obs.(exp.) 


Signal Region 

Background 

Observation 

c95 

^ obs.(exp.) 

(Jieo-) 

oL(exp.) [fb] 

SRHH 

3.1 ± 1.2 

3 

5.5 (5.5:;i) 

0.27 

(»47:S:^) 

SRLM 

22.1+4.7 

20 

12.4(13,2:«) 

0.61 

(0.65«;“) 

SRHM 

2.1 ± 1.5 

3 

6.4(5.2«;9 

0.31 



the expected background and the background uncertainty as input to the calculation. Possible signal 
contamination in the control regions is neglected. Dividing the limits on the number of BSM events by 
the integrated luminosity of the data sample, these can be interpreted as upper limits on the visible BSM 
cross section, cTyis = cr x Jix e, where cr is the production cross section for the BSM signal, Jl is the 
acceptance defined as the fraction of events passing the geometric and kinematic selections at particle 
level, and e is the detector reconstruction, identification and trigger efficiency. Table 8 summarises, for 
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Table 9: Acceptance times efficiency for the various signal regions, for a few selected (scalar top, scalar tau) signal 
mass hypotheses. For each mass point, values are shown only for channels targeting that point. The lepton-lepton 
results are taken from Ref. [50]. 


t\ mass 
[GeV] 

f 1 mass 
[GeV] 

lepton-lepton 

Axe 

lepton-hadron 
Axe (SREM) 

lepton-hadron 
Axe (SRHM) 

hadron-hadron 

Axe 

153 

87 

- 

1.29 X 10-4 

- 

2.27 X 10-4 

195 

87 

- 

1.36 X 10-4 

- 

4.46 X 10-4 

195 

148 

1.71 X 10-^ 

7.80 X 10-5 

- 

7.00 X 10-4 

195 

185 

8.01 X 10“^ 

- 

- 

- 

391 

148 

7.32 X 10-4 

- 

9.44 X 10-4 

3.40 X 10-5 

503 

493 

1.03 X 10-2 

- 

- 

- 

561 

87 

- 

- 

1.74 X 10-5 

6.70 X 10-5 

561 

337 

- 

- 

1.30 X 10-2 

9.90 X 10-5 

561 

500 

- 

- 

8.68 X 10-5 

2.50 X 10-5 


each SR, the estimated SM background yields, the observed numbers of events, and the expected and 
observed upper limits on event yields from a BSM signal and on cTyis. Table 9 summarises, for each SR, 
the acceptance times efficiency for the relevant final state under various signal mass hypotheses. 

Exclusion limits are derived for the scalar top pair production, assuming the t\ decays with 100% BR 
into bvrf\ and the fi decays into a tau lepton and a gravitino. The fit used for these limits is similar 
to that described in Sect. 5, but it now includes the expected signal in the likelihood, with an overall 
signal-strength parameter constrained to be positive. The CRs and SRs are fit simultaneously, taking 
into account the experimental and theoretical systematic uncertainties as nuisance parameters. The signal 
contamination in the CRs is also taken into account. Exclusion contours are set in the plane defined by 
the t\ and f i masses. 

Systematic uncertainties on the signal expectations stemming from detector effects are included in the fit 
in the same way as for the backgrounds. Systematic uncertainties on the signal cross section due to the 
choice of renormalisation and factorisation scales and PDE uncertainties are calculated as described in 
Sect. 6. Unlike other nuisance parameters, the signal cross-section uncertainties are only used to assess 
the impact of a ±lcr variation on the observed limit. 

Eor each mass hypothesis, the expected limits are calculated for the hadron-hadron selection, the two 
lepton-hadron selections, and the statistical combination of the lepton-lepton selections described in 
Ref. [50]. The selection giving the best expected sensitivity is used to compute the expected and observed 
CEi value. The resulting exclusion contours are shown in Eig. 9. The limits for each individual channel 
are reported in Eig. 10. The black dashed and red solid lines show the 95% CE expected and observed 
limits, respectively, including all uncertainties except for the theoretical signal cross-section uncertainty 
(PDE and scale). The yellow bands around the expected limits show the ±lcr expectations. The red dotted 
±lcr lines around the observed limit represent the results obtained when varying the nominal signal cross 
section up or down by its theoretical uncertainty. Numerical limits quoted on the particle masses are taken 
from these -Icr theoretical lines. 

As can be seen from Eig. 9, models with a scalar top mass below 490 GeV are excluded. Depending on 
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the scalar tau mass, some models with scalar top masses up to 650 GeV are also excluded. The scalar 
top masses below 150 GeV are not fully considered but they are unlikely to be viable because the cross 
section times branching ratio for tit\ bTbr+X is more than 25 times larger than the cross section times 
branching ratio for the production of tt decaying into the same di-tau final state, and measurements of the 
tt cross section in various final sfafes [107-110] are in good agreemenf wifh fhe SM predicfion. 


t,t,* production, BR(t, ^ T,bv) = 1, BR(x^^ x G) = 1 



[GeV] 


Figure 9: Observed and expected exclusion contours at 95% CL in the (fi,fi) mass plane from the combination of 
all selections. The dashed and solid lines show the 95% CL expected and observed limits, respectively, including 
all uncertainties except for the theoretical signal cross-section uncertainty (PDF and scale). The band around the 
expected limit shows the +lcr expectation. The dotted +lcr lines around the observed limit represent the results 
obtained when varying the nominal signal cross section up or down by the theoretical uncertainty. The LEP limit 
on the mass of the scalar tau is also shown. 
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8 Conclusion 


A search for direct pair production of supersymmetric partners of the top quark decaying via a scalar tau 
to a nearly massless gravitino has been performed using 20 fb“^ of pp collision data at = 8 TeV, 
collected by the ATLAS experiment at the LHC in 2012. Scalar top candidates are searched for in events 
with either two hadronically decaying taus, one hadronically decaying tau and one light lepton, or two 
light leptons. Good agreement is observed between the Standard Model background estimate and the 
data. The first results from a hadron collider search for the three-body decay mode to the scalar tau are 
presented. In the context of the model considered, lower limits on the scalar top mass are set at 95% 
confidence level, and found to be between 490 GeV and 650 GeV for scalar tau masses ranging from the 
LEP limit to the scalar top mass. 


t,t,* production, BR{t, -»tibv) = 1, BR{t -»t G) = 1 





700 

mj [GeV] 


700 

trij [GeV] 


Figure 10; Observed and expected exclusion contours at 95% CL in the (fi,fi) mass plane from the hadron- 
hadron (top left), the lepton-hadron low-mass (top right), the lepton-hadron high-mass (bottom left) and the lepton- 
lepton selections of Ref. [50] (bottom right). The dashed and solid lines show the 95% CL expected and observed 
limits, respectively, including all uncertainties except for the theoretical signal cross-section uncertainty (PDF and 
scale). The band around the expected limit shows the ±lcr expectation. The dotted ±lcr lines around the observed 
limit represent the results obtained when varying the nominal signal cross section up or down by the theoretical 
uncertainty. The LEP limit on the mass of the scalar tau is also shown. 
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